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Light hadron spectrum and quark masses 

T. Kaneko^ 

'^High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 

Recent developments in lattice QCD calculations of the light hadron spectrum and quark masses are reviewed. 



1. Introduction 

One of the main goals of lattice QCD is to con- 
firm the validity of QCD as the theory of strong 
interactions. Calculation of the hadron spectrum 
is fundamental in this respect. In the quenched 
approximation, large effort over the years has cul- 
minated in the work by CP-PACS|] which re- 
vealed a systematic deviation of the quenched 
spectrum from experiment. Since then, the fo- 
cus has shifted to full QCD including dynamical 
quarks, and a number of large-scale simulations 
have been executed in Nf ~ 2 full QCD|^ . A cru- 
cial issue here is how sea quark effects manifest 
in the spectrum. Work to explore this question 
has continued this year. 

A realistic QCD simulation has to treat the 
strange quark dynamically. The algorithmic de- 
velopment to treat odd number of flavors has ac- 
celerated in the last few years This has led to 
several attempts toward realistic calculations in 
N f — 3 QCD very recently. 

One of the major uncertainties in lattice cal- 
culations arises from the chiral extrapolation to- 
ward realistically light quark masses. One possi- 
bility to control this extrapolation is to use chiral 
perturbation theory (ChPT)|^ as a guide. Test 
of the validity of ChPT is an important step for 
this purpose, and full QCD data to examine this 
issue are becoming available. 

Another recent trend is the improvement of for- 
mulations of lattice QCD such as the domain- wall 
fcrmion|5| and fermions defined on anisotropic 
lattices [g. Quenched QCD presents a testing 
ground of these formulations, and some realistic 
simulations have been made. 

Spectrum calculations allow a simultaneous de- 
termination of the strong coupling constant and 



quark masses which are the fundamental param- 
eters of the Standard Model. Precise determina- 
tion of these quantities is an important issue in 
lattice QCD|Q,||. New studies in this area have 
also been reported at this conference. 

In this review, we present the status of lat- 
tice calculations of the light hadron spectrum and 
the fundamental parameters of QCD with focus 
on the points above. We concentrate on results 
obtained with the Wilson-type quark actions, as 
those with the Kogut-Susskind (KS) action are 
covered in a separate talkj|]. In Sec. ^, we dis- 
cuss sea quark effects in iV/ = 2 full QCD. Recent 
developments toward Nf — 3 simulations are re- 
viewed in Sec. Is], and verification of ChPT is dis- 
cussed in Sec. Section || is devoted to recent de- 
velopments of spectroscopic studies in quenched 
QCD. The status of the strong coupling constant 
and quark masses are updated in Sees. ^ and |^. 
A brief conclusion is given in Sec. ^j. 

2. Light hadron spectrum in Nf — 2 QCD 

2.1. Recent simulations 

Recent simulations in two-flavor QCD are listed 
in Table Most of the simulations with the 
Wilson-type quark action were made with the 
plaquette gauge action at a single lattice spac- 
ing a ~ 0.1 fm. An exception is the CP-PACS 
study l2l| which explored the range a~0.1-0.2 fm 
with the use of an improved gauge action. This 
work reported clear sea quark effects in the meson 
spectrum after the continuum extrapolation. 

New results were reported by UKQCD|i^ and 
JLQCD|i^ at this conference. The two simu- 
lations still use the plaquette gauge action and 
work at a single lattice spacing a ^ 0.1 fm, but 
the quark action is fully 0{a) improved with csw 



Table 1 

Recent simulations in Nf — 2 QCD. Values of csw for the clover quark action jio) are shown in brackets, 
where NP and TP stand for non-perturbative and tadpole-improved values, respectively. 



group 



gauge 



quark 



fls [fm] 



Ls[im\ mps,soa/™v,sca ref. 



SESAM 
TxL 

SESAM-TxL 



plaquette 
plaquette 
plaquette 



Wilson 
Wilson 
Wilson 



0.08 
0.08 
0.09 



1.3 
1.9 
1.5 



0.68-0.83 
0.57, 0.70 
0.68-0.85 



[^1 
[12 

[13 



UKQCD 

UKQCD 

UKQCD 

QCDSF 

JLQCD 



plaquette 
plaquette 
plaquette 
plaquette 
plaquette 



clover (1.76) 
clover (NP) 
clover (NP) 
clover (NP) 
clover (NP) 



0.12 
0.10 
0.10 
0.09 
0.09 



1.0- 1.9 
1.7 

1.6 

2.1,1.5 

1.1- 1.8 



0.67-0.86 
0.70-0.84 
0.58 

0.69, 0.76 
0.60-0.80 



CP-PACS RG-improvedQ clover (TP) 0.11-0.22 2.5-2.6 0.55-0.81 



[14 
[15 
[15 

[LI 

[L8jl|]_ 



[21 



Columbia 

MILC 

Columbia 



plaquette 
plaquette 
plaquette 



KS 
KS 

KS (g^ 1.8-5.0) 



0.09 

0.10-0.32 
0.23-0.34 



1.5 

2.4-3.8 
3.7-5.4 



0.57-0.70 

0.3-0.8 

0.35-0.80 



[22 
[23 
[24 



MILC 



Symanzik[25 



improved KS[ 26 1 0.13 



2.6 



0.50 



[27 



determined non-perturvatively [p8| . 

The UKQCD simulation shifts (3 with sea quark 
mass, keeping the Sommer scale ro/ap9[| fixed. 
The JLQCD runs were performed at a fixed 
f3. Another important difference is the range 
of quark mass covered in the two simulations: 
JLQCD explored light sea quark masses down to 
"^PS,sea/wv,sea~0.6, whilst the UKQCD's hght- 
est point is around mps.soa/^v.sca 0.7.[] We 
note that Orth et al. attempts to reduce quark 
masses to even lighter points, but actual simula- 
tions have not yet been started[^. 

2.2. Sea quark effects in meson spectrum 

Figure [l] shows how the valence quark mass 
dependence of the vector meson mass varies as 
a function of sea quark mass in the JLQCD 
and UKQCD simulations. In both UKQCD and 
JLQCD data, the set of points for the sea quark 
masses down to TOps.sca/™v,sca~0.7 almost over- 
lap. On the other hand, the JLQCD data for 
a lighter sea quark TOps,sca/'Tiv,soa ~ 0.6 clearly 
show a larger slope, and are closer to the experi- 
mental points marked by asterisks. 

The smallest sea quark mass in most of the 
previous simulations stopped at mps.sca/'riv.sca'^ 

t Although UKQCD made another simulation at a smaller 
sea quark mass (r7ips_aca/'"v,sca ~ 0.58, 1.6 fm), 

finite size effects seem to be significant there [hg||. 



0.7, which explains why sea quark effects were 
not clearly seen in these simulations. The 
CP-PACS data did cover the range down to 
™ps,sca/'7tv,sGa ~ 0.6. The trend of sea quark ef- 
fects as observed by JLQCD is present at a ~ 
0.1 fm, while it is less clear on coarser lattices 
due to larger error of ropl|. 

Figure y shows recent results of the J param- 
eter [^ calculated with differentiation in the va- 
lence quark mass as a function of the sea quark 
mass. The JLQCD data show a trend that 
the value of J, which is consistent with that of 
quenched QCD for heavier sea quark masses, in- 
creases as the sea quark mass decreases. The lat- 
ter feature can be observed also in the UKQCD 
data. As a result, the values extrapolated to the 
physical sea quark mass are closer to experiment 
than in quenched QCD. 

The J parameter in full QCD can also be 
obtained using the experimental definition J = 
mK'{rnK» —Tnp)/{rn\ — rn^), or from the com- 
bined chiral extrapolation as a function of sea 
and valence quark masses. The results are also 

plotted at (rgTOpS.sca 

)2 = in Fig. |. They lead 
to the same conclusion that J increases with dy- 
namical quarks. An exception is the SESAM 
point presumably due to heavier sea quark masses 
'7ips,sca/'7iv,sca^0.7 and scaling violation in their 
simulations. Recent results with KS fermions also 
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Figure 1. Vector meson mass as a function of 
pseudo-scalar meson mass squared, normalized in 
units of ro = 0.49 fm from JLQCD and UKQCD. 
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Figure 2. Recent results of J parameters in 
Nf = 2 full QCD and quenched QCD. Results 
presented this year are plotted in filled symbols. 

observed similar sea quark effects p7|. 

The results so far refer to a single lattice spac- 
ing of a ~ 0.1 fm. In Fig. |^ we reproduce the 
previous results by CP-PACSj2^ on the strange 
vector meson masses for a set of lattice spacings, 
and superimpose the new JLQCD points. The 
JLQCD points are consistent with the CP-PACS 
points both for Nf = 2 and Nf = (quenched), 
showing an increase of hyperfine splitting as dy- 
namical quark effects are included. 

Strictly speaking comparison of the two data 
should be performed in the continuum limit, al- 
beit scaling violation in the JLQCD data may 
be expected to be small by the use of the non- 
perturbatively 0(a)-improved action. This point 
should be checked in future studies. 



a [GeV"'] 

Figure 3. Strange vector meson masses with rriK 
as input in Nf — 2(filled symbols) and quenched 
QCD (open symbols). 

2.3. Baryons and finite size eff"ects 

Significant sea quark effects are also expected 
in the baryon spectrum, as it sizably deviates 
from experiment in quenched QCD|0|. However, 
no group has found clear evidence. A possible 
reason is finite size effects, which we expect to be 
more important for baryons than for mesons. 

Suppression of systematic errors to a few per- 
cent level is required for a convincing study of 
sea quark effects. In order to estimate the lat- 
tice size needed to achieve this accuracy, we fit 
hadron mass data of UKQCD and JLQCD (l9) 
on several lattice volumes to an ansats m{L) = 
m{L = oo) + c/i'^Q. Using the fitted value of 
m{L — oo), we show in Fig. ^ the relative magni- 
tude of finite size effects as a function of the spa- 
tial extent L (divided by tq) at sea quark masses 
corresponding to mps/my — 0.7 and 0.6. 

We observe that a lattice size of L ~ 2 fm suf- 
fices for mesons to keep finite size effects under 
a few percent level for sea quark masses down to 
'Ttps/wy ~ 0.6. On the other hand, calculations 
of the baryon masses need L ~ 3 fm at the same 
sea quark masses and the required size increases 
rapidly as the sea quark mass decreases. Further 
studies on such large volumes will be necessary to 
identify sea quark effects on the baryon spectrum. 

2.4. Glueball masses 

An unambiguous experimental identification of 
the glueball is still missing, and lattice QCD can 
provide helpful theoretical estimates. Calcula- 
tions in quenched QCD have been performed for 
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Figure 4. Mass shift at finite lattice size L as 
a function of L/tq. The dashed and dotted lines 
shows fit curve for data at sea quark mass corre- 
sponding to mps/TOv = 0.7 and 0.6, respectively. 

a long time and recent results for J^^ — 0++ 
and 2++ states are in good agreement with each 
other iQ. The focus has now shifted to calcula- 
tions in full QCD and mixings with quarkonia. 

UKQCD|| and SESAM-TxL|3t| reported 
new results this year. Both groups find a decrease 
of mo++ by about 20% with dynamical quarks. 
The observed sea quark mass dependence, how- 
ever, is quite different: SESAM-T^L finds a sig- 
nificant slope, while UKQCD data are almost flat. 
It is, thus, too early to draw a definite conclusion 
on the sea quark effect. Scaling violation should 
be also studied carefully, as it is known to be large 
in quenched QCD. A similar indication was al- 
ready seen by SESAM-TxL[|6) in full QCD. 

3. QCD with dynamical up, down and 
strange quarks 

Realistic simulations of QCD require up, down 
and strange quarks to be treated dynamically. 
While such studies already exist for the KS quark 
action using an approximate R algorithmic,^ , 
exact algorithms for odd number of Wilson-type 
quarks were developed only recently 

This year JLQCD initiated a study of Nf ~ 3 
full QCD using an exact algorithm for the 0(a)- 
improved Wilson fermion|^ . As a first step, they 
studied the phase structure on the (/3, k) plane 
and found an unexpected first-order phase tran- 
sition with the standard plaquette gauge and the 
0(a)-improved Wilson quark action, as shown in 
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Figure 5. Phase diagram in Nf — 3 full QCD 
with the plaquette gauge and 0(a)-improved Wil- 
son quark. 

Fig. Hl^^ . The transition points do not move be- 
tween 8"^ X 16 and 12^ x 32 lattices, and they dis- 
appear at /3>5.0. Hence this is a bulk transition 
at zero temperature, and is a lattice artifact. 

Their finding poses a serious practical prob- 
lem, because the lattice spacing estimated from 
To equals ^ 2.6 GeV at the end-point of the 
first-order transition at /3~5.0. Large-scale simu- 
lations at such fine lattice spacings are not practi- 
cal with the computer power currently available. 

They found, however, that this unphysical 
phase transition disappears if improved gauge ac- 
tions are employed. For realistic simulations, it 
is a natural choice to use improved gauge actions 
at moderately large values of lattice spacing such 
as a ~ 0.1-0.2 fm. 

4. Test of chiral perturbation theory 

Chiral perturbation theory (ChPT)0| predicts 
the presence of logarithm singularities in the 
quark mass dependence of physical quantities, 
e.g., for the pseudo-scalar meson mass mpg, 
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where y = 2Boniq/{ATrf)^. A test of this relation 
is shown in Fig. |g, where lattice data for two- 
flavor QCD from JLQCD ||l9| are plotted together 
with the fit curves of Eq. ^ ) assuming / to be a 
free parameter or fixed to / = 93 MeV. The lat- 
tice results are not consistent with the curvature 
characterized by the chiral logarithm. 
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Figure 6. Test of chiral perturbation theory. 

Recently UKQCD|39j attempted to extract the 
low energy constants of chiral lagrangian from a 
fit of lattice data to the partially quenched ChPT 
formula pof using ratios proposed by ALPHA p|. 
The curvature expected from the chiral logarithm 
is also not clear in their data. 

A possible reason behind these disagreements 
is that sea quark masses in current simulations 
are too large. In fact, the coefficient of the chiral 
logarithm is modified when the effect of flavor sin- 
glet meson is included into ChPT[Q; it can take 
much smaller values than 1/Nf, if the pesudo- 
scalar meson is not so light compared to the fla- 
vor singlet meson ll^. A recent work by Nelson et 
a/.p3|] using KS fermions suggests a similar con- 
clusion. The curvature is visible in their data at 
small sea quark masses, which are possible with 
KS fermions. 

Another interesting subject of ChPT is deter- 
mination of the low energy constants ai from lat- 
tice data |^ , ^ . In particular, ag is an important 
quantity since it is related to the possibility of 
TO„ = 0j4j]. UKQCD attempted to extract a^'s in 
two-flavor QCD, and obtained as = 0.79(17) (21), 
which disfavors to„ = 0. Nelson et al. obtained a 
consistent conclusion in three flavor QCD. These 
results, however, should be confirmed with lighter 
sea quark masses, for which ChPT at the one-loop 
level becomes more reliable. 

5. Developments in quenched QCD 

5.1. Test of improved formulations 

Quenched QCD is now used to develop and test 
improved formulations of lattice QCD. Domain- 
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Figure 7. Continuum extrapolation of m/^ . /rup 
obtained with domain-wall fermion(squares) and 
on anisotropic lattice(triangles). 

wall fermion formalism is an important case, as 
it is expected to realize exact chiral symmetry at 
finite lattice spacings. Anisotropic lattices may 
provide a way to achieve simulations with large 
lattice cut-offs on large spatial volumes. 

Figure shows results of m^c . /nip mass ratio 
calculated with the domain- wall fermions ||4^, and 
on anisotropic lattices Clearly these formu- 
lations reproduce, and confirm, the small hyper- 
fine splitting in the quenched meson spectrum ob- 
tained by the precison calculation by CP-PACS 
using the conventional Wilson quark action[0. 

Recently, the Wilson quark action with a chi- 
rally twisted mass term psf has drawn attention 
as an alternative formulation to suppress excep- 
tional configurations . Some evidences of such 
suppression have been presented at this confer- 
ence^,^. It will be interesting to apply this 
action to the non-perterburtive determination of 
improvement and renormalization constants at 
strong coupling, where exceptional configurations 
cause a serious problem i 
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5.2. Negative parity baryon 

An interesting subject in recent quenched sim- 
ulations is to explore the rich structure of the ex- 
cited baryons. There already exist several stud- 
ies for the negative parity baryon N* [J^ = 
1/2^) at finite lattice spacings fs^. The first 
systematic study of scaling violation and finite 
size effects has been performed this year by the 
QCDSF-UKQCD-LHPC collaboration using the 
0(a)-improved Wilson action[E5|. Their result 



totv* /"T-Af = 1-50(3)(8)(15) in the continuum limit 
is consistent with the experimental value 1.63. 

It is also interesting to explore other excita- 
tions. A preliminary result for the excited decu- 
plet baryons, J^ = 3/2+ and 3/2 ~, was presented 
by Lee et al. at this conference!^. 

6. Strong coupling constant 

Recent results of a"J^~^\Mz) from lattice 

MS ^ ^ 



Table 2 

Recent results of ai^~^^ {Mz)- AH results except 
Ref. [pO| are calculated through ap [M . 



QCD are summarized in Table Recent val- 
ues using the Wilson-type quark action for sea 
quarks | |58| , |6^ are significantly smaller compared 
to a standard lattice result previously obtained 
by Davies et aL[|6l| using KS fermions. In par- 
ticular this year's result of the QCDSF-UKQCD 
collaboration is smaller despite a consistent use 
of two-loop perturbation theory and a continuum 
extrapolation . 

The origin of the difference between results 
with the KS and Wilson fermions should be clar- 
ified. Since the discrepancy is already present in 
the value of ap at the cutoff scale in two- flavor 
QCD, systematic uncertainties such as scaling vi- 
olation and higher order corrections in perturba- 
tion theory should be studied carefully in future. 

Calculations of as through various definitions 
are also helpful to resolve this problem. Boucaud 
et al. presented a preliminary result from lat- 
tice calculations of the gluon green functions j60| . 
A non-perturbative determination of Ug in the 
Schrodinger functional scheme was performed 
by ALPHA[||| in two-fiavor QCD, and determi- 
nation of the lattice energy scale is in progress . 

7. Quark masses 

By last year results for the light quark mass 
mud in quenched QCD converged within an ac- 
curacy of about 15% [H. For nig there turned out 
an additional uncertainty of about 20% due to 
the choice of the input meson mass {e.g. vtik or 
m0)||l|. A noticable trend in this year's calcu- 
lation is the use of the domain-wall and overlap 
fermions. Calculations of the charm quark mass 
rric have also been pursued in quenched QCD. 

In iV/ = 2 fuU QCD, CP-PACS reported results 
of a systematic study already last year pa . They 



group 


quark al^ 


SESAM|8| 
QCDSF-UKQCD 1^ 
Boucaud et a/.pO|] 


Wilson 0.1118(17) 
1 clover(NP) 0.1076(20) (18) 
Wilson 0.113(3)(4) 


Davies et a/.[31| 


KS 0.1174(24) 



found that sea quark effects reduce rriud and rUg 
by about 25%, and the discrepancy in rus with 
vriK or as input is remarkably reduced. This 
year there has been less progress in full QCD. 

7.1. Strange quark mass in quenched QCD 

A new estimate of nig was obtained by CP- 
PACS using the domain- wall fermion|6^. Their 
calculations were performed at ^ 2 and 

3 GeV using one-loop perturbative renormal- 
ization constants. Giusti et al. calculated 
{nig -\- iriud) at a^^ ^ 2 GcV using the over- 
lap fermion]6^. Their rcnomalization constants 
are evaluated non-perturbatively in the RI/MOM 
scheme. 

Recent results with ttik as input are plotted in 
Fig. ||. A reasonable agreement among the new 
and previous results is observed. Therefore, the 
average mf^{2 GeV) = 105 ± 15 MeV does not 
change significantly since last yearj^. 

7.2. Light quark mass in quenched QCD 

We consider the ratio uig/niud, rather than 
niud itself, since systematic uncertainties partially 
cancel in this ratio. A new result mg/rriud = 
26.2(2.3) by CP-PACS[|6) as well as previous re- 
sults are in good agreement with the prediction 
of one-loop ChPT : 24.4(1. 5)|6|. We therefore 
apply the ratio from ChPT to the average of nig 
to obtain m^{2 GeV) = 4.3 ± 0.7 MeV. 

Dong et al. calculated niud with the over- 
lap fermion[p9[. Their estimate (2 GeV) = 
5.3(0.6) (0.7) MeV agrees with the above average. 

7.3. Charm quark mass in quenched QCD 

There have been few studies for calculating 
the charm quark mass TOc||7C|]. The reason is 
that 0((aTOq)") scaling violation is expected to 
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Figure 8. Continuum extrapolation of strange 
quark mass with rriK as input in MS scheme 
at scale /i = 2 GeV. Filled and open symbols 
show new and previous results, respectively. The 
choice of the quark action is written in brackets, 
where "W"=Wilson, "C"=clover, "KS"=Kogut- 
Susskind and "DWF" =domain-wall fermions. 

be large in relativistic formulations, while the 
charm quark is too light to be treated in the non- 
relativistic approximation. 

New estimates of rric have been reported by two 
groups this year, both using the 0(a)-improved 
action to reduce scaling violation. Becirevic et al. 
performed a simulation at a fixed lattice spacing 
~3 GeV0|. The results significantly differ 
for rric obtained from vector(VWI) and axial vec- 
tor Ward identities (AWI). They take an average 
to quote mf^{mc) = 1.26(3)(12). The difference 
between the two methods(AWI and VWI) and the 
average is treated as a systematic error. 

ALPHA performed simulations in a wide range 
of lattice spacing a^^ ^ 2-4 GeV and performed 
the continuum extrapolation |f72| . They observed 
a nice convergence of rric determined from AWI 
and VWI toward the continuum limit, and ob- 
tained m^(mc) = 1.314(40)(20)(7) in the con- 
tinuum limit. 

These results are plotted in Fig. ^ together with 
previous results |f70| and an average by PDG[f73|. 
The two new results are consistent with previous 
ones as well as the PDG's average. 

A new estimate using the heavy quark action 
in the FNAL interpretation |74| was presented by 
Juge et al. at this conference |75|]. While the esti- 
mation of the systematic error is still in progress. 
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Figure 9. Recent results of charm quark mass. 
An average quoted by PDGj?^ is also plotted. 

their preliminary result is also around 1.3 GeV 
and consistent with the above new results. 

7.4. Results in iV/ = 2 full QCD 

SESAM-T^L attempted a continuum extrapo- 
lation of rriud and nig from previous and new 
simulations at /3 = 5.5[^. Their range of the lat- 
tice spacing is, however, narrow and far from the 
continuum limit. This leads to a large uncertainty 
in their results, mjf|(2 GeV) =4.5(1.7) MeV and 
mf^{2 GeV) = 92(83) MeV. In near future, we 
expect results from the UKQCD and JLQCD sim- 
ulations with the 0(a)-improved Wilson fcrmion. 

8. Conclusions 

In quenched QCD, improved formulations, such 
as the domain-wall and overlap fermions and 
anisotropic lattice actions, have been tested by 
calculating the light hadron spectrum and quark 
masses, and good consistency with previous re- 
sults has been observed. These formulations are 
ready for applications to precise determinations 
of other observables. 

In Nf = 2 full QCD, sea quark effects on the me- 
son spectrum, which were previously observed by 
CP-PACS, have been confirmed using the 0(a)- 
improved Wilson and an improved KS fermions. 
However we still need further studies for sea quark 
effects on the baryon and glueball spectrum. Ex- 
ploring small sea quark masses with Wilson-type 
fermions is also an important issue to resolve the 
poor consistency between lattice data and ChPT. 

Pursuit of realistic simulations in three flavor 



QCD is one of the most important issue in spec- 
troscopic studies. This year significant progress 
to this end has been made by several groups us- 
ing both of the KS and Wilson fermions. Further 
systematic studies will be expected in near future. 

And finally, there has been few progress in non- 
perturbative calculation of the improvement co- 
efficient ca and renormalization factors in full 
QCD. These calculations arc urgently required for 
precise determinations of the quark masses and 
decay constants. 

Acknowledgments 

I thank Y. Aoki, G.S. Bah, C. Bernard, 
S. Collins, G. Di Carlo, S.J. Dong, N. Eicker, 
R. Frezzotti, C. Gcbcrt, L. Giusti, A. Hart, 
S. Hashimoto, S. Hauswirth, J. Hein, K. Ide, 
A.C. Irving, K-I. Ishikawa, K. Jansen, F.X. Lee, 
L. Lcvkova, X.Q. Luo, T. Manke, H. Matsu- 
furu, C. McNcilc, H. Moutardc, A. Nakamura, 
M. Okawa, K. Orginos, O. Philipscn, C. Rcbbi, 
J. Rolf, G. Rossi, S. Sint, H. Wittig, T. Yamazaki, 
and A. Ukawa for communications and useful dis- 
cussions. I also thank S. Hashimoto, Y. Iwasaki, 
M. Okawa and A. Ukawa for valuable suggestions 
on the manuscript. This work is supported by the 
JSPS Research Fellowship. 

REFERENCES 

1. S. Aoki et al. (CP-PACS Collaboration), 
Phys. Rev. Lett. 84 (2000) 238. 

2. For a recent review, see S. Aoki, Nucl. Phys. 
B (Proc. Suppl.) 94 (2001) 3. 

3. M. Peardon, review in these proceedings. 

4. J. Gasser and H. Leutwyler, Ann. Phys. 158 
(1984) 142; Nucl. Phys. B250 (1985) 465. 

5. Y. Kikukawa, review in those proceedings and 
references therein; P. Hernandez, ibid. 

6. T.R. Klassen, Nucl. Phys. B (Proc. Suppl.) 
73 (1999) 918. 

7. For a recent review, see J. Shigemitsu, Nucl. 
Phys. B (Proc. Suppl.) 53 (1997) 16. 

8. For a recent review, sec V. Lubicz, Nucl. 
Phys. B (Proc. Suppl.) 94 (2001) 116. 

9. D. Toussaint, review in these proceedings. 

10. B. Sheikholeslami and R. Wohlert, Nucl. 
Phys. B259 (1985) 572. 



11. N. Eicker et al. (SESAM Collaboration), 
Phys. Rev. D59 (1999) 014509. 

12. T. Lippert et al. (SESAM-TxL Collabora- 
tion), Nucl. Phys. B (Proc. Suppl.) 60A 
(1998) 311. 

13. N. Eicker (SESAM-TxL Collaboration), 
these proceedings. 

14. C.R. Allton et al. (UKQCD Collaboration), 
Phys. Rev. D60 (1999) 034507. 

15. A.C. Irving (UKQCD Collaboration), Nucl. 
Phys. B (Proc. Suppl.) 94 (2001) 242. 

16. C.R. Allton et al. (UKQCD Collaboration), 
hep-lat/0107021; these proceedings. 

17. D. Pleiter (QCDSF-UKQCD Collaboration), 
Nucl. Phys. B (Proc. Suppl.) 94 (2001) 265; 
H. Stiiben (QCDSF-UKQCD Collaboration), 
ihid. 94 (2001) 273. 

18. S. Aoki et al. (JLQCD Collaboration), Nucl. 
Phys. B (Proc. Suppl.) 94 (2001) 233. 

19. S. Hashimoto (JLQCD Collaboration), these 
proceedings. 

20. Y. Iwasaki, Nucl. Phys. B258 (1985) 141; 
Univ. of Tsukuba report UTHEP-118(1983), 
unpublished. 

21. A. Ah Khan et al. (CP-PACS Collaboration), 
Nucl. Phys. B (Proc. Suppl.) 94 (2001) 229; 
hep-lat/0105015. 

22. C. Sui et al, Nucl. Phys. B (Proc. Suppl.) 
73 (1999) 228; see also R.D. Mawhinney, ibid. 
83-84 (2000) 57. 

23. C. Bernard et al. (MILC Collaboration), 
Nucl. Phys. B (Proc. Suppl.) 73 (1999) 198; 
ibid. 60A (1998) 297 and references therein. 

24. L. Levkova, these proceedings. 

25. P. Weisz, Nucl. Phys. B212 (1983) 1; 
M. Liischer and P. Weisz, Comm. Math. Phys. 
97 (1985) 59. 

26. S. Naik, Nucl. Phys. B316 (1989) 238; 
K. Orginos, D. Toussaint and R.L. Sugar, 
Phys. Rev. D60 (1999) 054503; Nucl. Phys. 
(Proc. Suppl.) 83-84 (2000) 878. 

27. C. Bernard et al. (MILC Collaboration), 
Phys. Rev. D64 (2001) 054506. 

28. K. Jansen and R. Sommer (ALPHA Collabo- 
ration), Nucl. Phys. B530 (1998) 185. 

29. R. Sommer, Nucl. Phys. B411 (1994) 839. 

30. B. Orth, these proceedings. 

31. CP-PACS Collaboration, unpublished. 



32. P. Lacock and C. Michael, Phys. Rev. D52 
(1995) 5213. 

33. M. Fukugita et al, Phys. Lett. B294 (1992) 
380; Phys. Rev. D47 (1993) 4739. 

34. For a recent review, see D. Toussaint, Nucl. 
Phys. B (Proc. Suppl.) 83-84 (2000) 151. 

35. A. Hart and M. Teper (UKQCD Collabora- 
tion), hep-lat/0108022; A. Hart (UKQCD- 
QCDSF Collaboration), these proceedings. 

36. G.S. Bah et al. (SESAM-TxL CoUaboration), 
Phys. Rev. D62 (2000) 054503; G.S. Bah, hep- 
ph/0110254. 

37. K-1. Ishikawa (JLQCD Collaboration), these 
proceedings. 

38. M. Okawa (JLQCD CoUaboration), these 
proceedings. 

39. A.C. Irving et al. (UKQCD Collaboration), 
Phys. Lett. B518 (2001) 243. 

40. S.R. Sharpe, Phys. Rev. D56 (1997) 7052. 

41. J. Heitger, R. Sommer and H. Wittig (AL- 
PHA Collaboration), Nucl. Phys. B588 (2000) 
377. 

42. M.F.L. Goltcrniann and K.C. Leung, Phys. 
Rev. D57 (1998) 5703. 

43. D. Nelson, these proceedings. 

44. A.G. Cohen, D.B.Kaplan and A.E.Nelson, 
JHEP 11 (1999) 027. 

45. S. Sharpe and N. Shoresh, Nucl. Phys. B 
(Proc. Suppl.) 83-84 (2000) 968. 

46. Y. Aoki (RBC Collaboration), these proceed- 
ings. 

47. H. Matsufuru, T. Onogi and T. Umeda, hep- 

lat/0107001; these proceedings. 

48. S. Aoki, Phys. Rev. D30 (1984) 2653. 

49. W. Bardeen et al, Phys. Rev. D59 (1999) 
014507; R. Frezzotti et al, Nucl. Phys. B 
(Proc. Suppl.) 83-84 (2000) 941; JHEP 0108 
(2001) 058. 

50. R. Frezzotti et al, these proceedings. 

51. C. McNeile (UKQCD Collaboration), these 
proceedings. 

52. M. Liischer et al, Nucl. Phys. B491 (1997) 
323. 

53. K. Ide (CP-PACS Collaboration), these pro- 
ceedings. 

54. F.X. Lee and D.B. Lcinwcbcr, Nucl. Phys. 
B (Proc. Suppl.) 73 (1999) 258; F.X. Leo, 
ibid. 94 (2001) 251; S. Sasaki, T. Blum and 



S. Ohta, hep-lat/0102010. 

55. M. Gockeler et al (QCDSF-UKQCD-LHPC 
Collaboration) , hcp-lat /0106022. 

56. F.X. Lee, these proceedings. 

57. G.P. Lepage and P.B. Mackenzie, Phys. Rev. 
D48 (1993) 2250. 

58. A. Spitz et al (SESAM CoUaboration), Phys. 
Rev. D60 (1999) 074502. 

59. S. Booth et al (QCDSF-UKQCD CoUabora- 
tion), Phys. Lett. B519 (2001) 229; G. Schier- 
holz (QCDSF-UKQCD Collaboration), these 
proceedings. 

60. Ph. Boucaud et al, hep-ph/0107278; 
H. Moutarde, these proceedings. 

61. C.T.H. Davies et al, Phys. Rev. D56 (1997) 
2755. 

62. L. Marcantonio et al (UKQCD Collabora- 
tion), Nucl. Phys. B (Proc. Suppl.) 94 (2001) 
363. 

63. M. Liischer et al, Nucl. Phys. B413 (1994) 
481. 

64. A. Bode et al (ALPHA CoUaboration), Phys. 
Lett. B515 (2001) 49; U. Wolff (ALPHA Col- 
laboration), these proceedings. 

65. A. Ah Khan et al (CP-PACS CoUaboration), 
Phys. Rev. Lett. 85 (2000) 4674. 

66. A. Ah Khan et al (CP-PACS Collaboration), 
hep-lat/0105020. 

67. L. Giusti, C. Hoelbling and C. Rebbi, hep- 
lat/0108007; these proceedings. 

68. H. Leutwyler, Phys. Lett. B378 (1996) 313. 

69. S.J. Dong, these proceedings. 

70. A. Bochkarev and P. de Forcrand, Nucl. Phys. 
B477 (1996) 489; V.Gimenez et al, Nucl. 
Phys. B540 (1999) 472; A.S. Kronfeld, Nucl. 
Phys. B (Proc. Suppl.) 63 (1998) 311. 

71. D. Becirevic, V. Lubicz and G. Martinelli, 
hep-ph/0107124. 

72. J. Rolf and S. Sint (ALPHA Collaboration), 
these proceedings. 

73. D.E. Groom et al (Particle Data Group), 
Eur. Phys. J. C15 (2000) 1. 

74. A.X. El-Khadra, A.S. Kronfeld and 
P.B. Mackenzie, Phys. Rev. D 55 (1997) 3933. 

75. K.J. Jugc, these proceedings. 

76. N. Eicker et al (SESAM CoUaboration), 
Phys. Lett. B407 (1997) 290. 



